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Abstract. The local velocity field of young stars is domi- 
nated by the galactic rotation, the kinematics of the Gould 
Belt and the nearest OB associations and open clusters, 
and the kinematics of the spiral structure. We re-examined 
here this local velocity field by using a large sample of 
nearby O and B stars from the Hipparcos Catalogue. The 
high quality astrometric data are complemented with a 
careful compilation of radial velocities and Stromgren pho- 
tometry, which allows individual photometric distances 
and ages to be derived. The Gould Belt extends up to 
600 pc from the Sun with an inclination with respect to 
the galactic plane of = 16-22° and the ascending node 
placed at f^o = 275-295°. Approximately 60% of the stars 
younger than 60 Myr belong to this structure. The val- 
ues found for the Oort constants when different samples 
selected by age or distance were used allowed us to inter- 
pret the systematic trends observed as signatures induced 
by the kinematic behaviour of the Gould Belt. The con- 
tribution of Sco-Cen and Ori OBI complexes in the char- 
acterization of the expansion of the Gould Belt system is 
also discussed. We found that a positive K-teim remains 
when these aggregates are excluded. From the kinematic 
behaviour of the stars and their spatial distribution we 
derive an age for the Gould Belt system in the interval 
30-60 Myr. 

Key words: Galaxy: kinematics and dynamics - Galaxy: 
solar neighbourhood - Galaxy: structure - Stars: early- 
type - Stars: kinematics 



1. Introduction 

The kinematic study of the local system of young stars of- 
fers an excellent opportunity for understanding the history 
of recent star formation and for improving our knowledge 
of the dynamics involved in the evolution of our galaxy. 
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The motion of young nearby stars deviates consider- 
ably from the general field of galactic rotation (see e.g. 

1973j ). The presence of a positive 



du Mont 1977: Clube 



K-term, corresponding to an overall expansion of the lo- 
cal system of the earliest stars was long ago (Campbell 



1913) recognized as the main kinematic characteristic of 



the Gould Belt. This system is recognized to be a great 
circle inclined some 20° respect to the galactic equator 
and is traced by young stars and OB associations, HI, 
molecular clouds and dust. A detailed review of the struc- 
ture, kinematics and origin of the G ould Be lt has recently 
been undertaken by Poppel ( 1997 ). Lesh ( 1968 ), Westin 



( 1985 ) and Comeron et al. ( 1994 ; hereafter referred to as 
CTG) found, in addition to the positive K-valae, depar- 
tures from the values of other Oort constants when work- 
ing with stars associated with the Gould Belt. In addition 
to the peculiar kinematics of the Gould Belt system, de- 
viations from circular motion have also been analysed in 
the context of the kin emati c effects of the s piral structure 
(Cre ze fc Menness ier |1973| ; Lindblad |l980t Byl & Oven- 



den 



1981; Westin 1985; Comeron & Torra 1991 



among 



others) 



The publication of the Hipparcos data (ESA 1997| ) has 
made it possible to re-analyse the galactic vel ocity field 
and its local irregularities. Feast & Whitelock ( |1997 ) re- 
determined the mean values of the Oort constants from 
Hipparcos proper motion of Cepheids (A = 14.8 ± 0.8 km 
s " 1 kp c" 1 , B = -12.4 ± 0.6 km s" 1 kpc" 1 ). Feast et al. 
( |1998| ) obtained A = 15.1 ± 0.3 km s" 1 kpc" 1 by using 
radial velocities and the new zero-point period-luminosity 
relation from Hipparcos trigonometric parallaxes. Lind- 



blad et al. (1997) and Torra et al. (1997) presented the 
first results concerning the structure and kinematics of 
the local system of young stars. Whereas Lindblad et al. 



(1997) do not rule out the possibility that the Gould Belt 
is just a random configuration of two or three dominat- 



ing associations, Sterzik et al. (199?) and Guillout et al. 
(19981), from a RASS-Tycho sample, suggest that the 
Belt is a disk-like rather than a ring-like structure. The 
Hippa rcos census of nearby OB associations (de Zeeuw 
et al. 1999 ) and the review of the mean astrometric pa- 
rameters of open clusters with Hipparcos (Robichon et al. 
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1999) have substantially increased our knowledge of the 
kinematic behaviour of the young stellar system in the 
solar neighbourhood. 

In Sect. |]we describe the samples we have used, based 
on the astrometric Hipparcos data plus a careful compi- 
lation of available radial velocities and Stromgren pho- 
tometry. This information enables us to obtain reliable 
space velocities and individual ages for a large number 
of stars. In Sect. |^ a maximum likelihood method is ap- 
plied in deriving the structural parameters and age of the 
Gould Belt. The velocity field was studied by means of 
the classical first-order approach in Sect. |J. The trends 
observed in the Oort constants are interpreted in terms 
of the expansion of the Gould Belt and the influence of 
stellar aggregates. 

2. The working sample 

Our initial sample (see Fernandez [1998 for more details) 
contained 6922 O- and B-type stars (Hipparcos Internal 
Proposal INCA060 completed with all the O and B survey 
stars), of which 5846 belong to the Hipparcos survey. The 
observational data were taken from the following sources: 

— Astrometric data from the Hipparcos Catalogue (ESA 



1997): positions in equatorial coordinates, parallaxes 



and proper motions, together with their standard er- 
rors and the correlations between them. Given the 
mean standard error in the trigonometric parallax pro- 
vided by Hipparcos, reliable distances are available 
only up to 200-400 pc. 

— Stromgren photometry from Hauck & Mermilliod's 
( |1998| ) compilation for deriving individual photometric 
distances and ages. 

— Radial velocities from Grenier's ( 1997| ) compilation 
plus additional sources. 

In the following subsections we describe the procedure 
used and the accuracy achieved in the derivation of indi- 
vidual distances, spatial velocities and ages, together with 
a discussion of the possible observational biases present in 
the final working samples. 

2.1. Stellar distances 



To derive the best distance estimate for each star in the 
sample, an analysis of individual errors, possible observa- 
tional biases and systematic differences between distances 
derived from Hipparcos trigonometric parallaxes and pho- 
tometric absolute magnitudes was performed. 

Only 3031 stars from our initial sample have complete 
Stromgren photometry (b — y, mi, c±, (3 and V) in Hauck 
& Mermilliod's (1998) catalogue. Recently, Kaltcheva & 
Knude ( 1998Q found that photometric distances derived 
from Crawford's (|1978|) and Balona & Shobbrook's ([19840 
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calibrations show good agreement with Hipparcos trigono- 
metric distances, and that no dependence on star rotation 



Fig. 1. Comparison of photometric and trigonometric dis- 
tances for the sample stars with o-^/tt < 0.15. The solid 
line is a 25-point running average. 



can be observed. The comparison of distances derived us- 
ing Crawford's calibration and Hipparcos distances for the 
non- binary stars in our sample with u^/tt < 0.15 (see Fig. 
[l]) shows no systematic trends. We also verify that the use 
of another calibration (Balona & Shobbrook 1984; Jakob- 
sen |1985 ) does not alter the kinematic results presented in 
Sect. After this analysis, Crawford's (1978) calibration 
was adopted to derive photometric distances. A relative 
error in the p hotom etric distance was computed follow- 
ing Lindroos ( |l98l| ). Depending on the spectral type and 
luminosity class, this error ranges between 14-23%. 

Reliable photometric distances can only be derived for 
single, non-variable and non-peculiar stars, so no photo- 
metric distance was derived for those stars classified as 
double or multiple in the Hipparcos Catalogue with a com- 
ponent separation p < 10" and a magnitude difference 
between components A.H p < 3 m , for variable stars with a 
variation in the Hipparcos magnitude system AH P > 0.6 m 
and for stars photometrically classified as peculiar (Jordi 
et al. |1997| ). 

For those stars for which we could obtain trigonometric 
and photometric distances, we used the distance with the 
smallest relative error. This procedure was preferred to 
the derivation of an error-weighted mean of both distance 
measurements since the latter would only systematically 
reduce the error for those stars around 150-250 pc (as 
will be seen in Sect. |4|, this error is used as a weight in 
the condition equations, so a different weight would be 
assigned as a function of distance). Furthermore, negative 
trigonometric parallax or the bias in the trigonometric 
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distance discussed below - non negligible for stars with 
large relative error in parallax - could not be controlled. 

In the case of stars for which the photometric distance 
was not available, to avoid bias, the trigonometric distance 
was only accepted if its relative error was smaller than 
25%. If distances are estimated as R = 1/tt, a symmetric 
error law for parallaxes results in a non-symmetric, biased 
distribution for distances. According to Arenou & Luri 
( |1999D , for small relative errors (Sa 25%) and assuming a 
Gaussian law for the error on the observed parallax, this 
bias can be approximated by: 



B(R) 



it \n t 



(1) 



being 7Tt the true parallax. Individual corrections for the 
observed parallaxes are not possible because the bias is a 
function of the true parallax (Brown et al. |1997| ). However, 
an estimation of the effect of this bias in our trigonometric 
distances may be made considering in Eq. (|l|) the observed 
instead of the true parallaxes. In Fig. || we show the rel- 
ative bias distribution (B(R)/R) for the 858 stars in our 
sample for which we chose the trigonometric distance as 
the best distance estimate. This relative bias is always 
less than 5.5%, and is smaller than 3% for 88% of the 
stars, thus giving a bias smaller than 5 pc for 82% of the 
stars. Therefore, given the impossibility of an individual 
correction of the biases, while recognizing their smallness 
(comparing, for instance, with the relative errors in the 
parallax), we may conclude that our assumption of con- 
sidering directly the trigonometric distance as given by 
R = 1/n for stars with relative error in parallax smaller 
than 25% is a good approximation. 



Fig. 3. Fraction of stars with radial velocity (top) and 
distribution of the stars with proper motion (blank his- 
togram) and radial velocity (filled histogram) as a func- 
tion of proper motion (bottom). Error bars were estimated 
from a Poissonian error distribution. 



2.2. Stellar radial velocities 

Our main source of radial velocities was the compilation 



of Grenier (1997), who made a complete compendium and 
revision of the compilations of Barbier-Brossat (1997) and 
Duflot et al. (1995). Priority was given to Barbier-Brossat, 
and only stars with A, B or C quality in Duflot et al. were 
considered. Using these sources, 3397 stars from our initial 
sample have radial velocity measurements. We rejected 
those stars with an individual error in the radial velocity 
higher than 10 km s _1 (131 stars, i.e. 3.9% of the stars 
with radial velocity). 



Binney et al. (1998), when working with nearby stars 
from all spectral types, emphasized that, due to obser- 
vational programmes, radial velocity availability is higher 
for high-proper motion stars. In our case, specific obser- 
vational programmes were undertaken in parallel with the 
Hipparcos mission to obtain radial velocity data for sur- 
vey early type Hipparcos stars. To evaluate the effects of 
the observational constraints on our kinematical study, we 
plotted in Fig. ||the fraction of stars with known radial ve- 
locities (qv T ) against the total proper motion. We can see 
as this fraction is not a flat function of it decreases 
for Ss 10 mas yr _1 and rises for 10 mas yr -1 . 

To understand this effect in Fig. [I| qv r is shown against 
proper motion and distance. We can see a higher degree 
of completeness for distant stars. From Fig. |^ we derived 
a completeness limit of V w 6.5 for radial velocities, al- 
though nearly all far and faint stars (R 1000 pc, so low 
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fi [mas/yr] 

Fig. 4. Fraction of stars with radial velocity against 
proper motion and distance. This fraction is showed in 
a grey scale, from (white) to 1 (black). 




Fig. 5. Fraction of stars with radial velocity against dis- 
tance and apparent visual magnitude. This fraction is 
showed in a grey scale, from (white) to 1 (black). 



have radial velocity data. A possible origin for this 
effect is due to observing programmes devoted to open 
clusters and associations. From this analysis we conclude 
that a kinematical bias is present in our sample and, al- 
though its effects on the analysis performed in Sect. [| are 
expected to b e negligible, they must be evaluated - see 
Appendix B: - through numerical simulations. 



2.3. Stellar ages 

Individual ages were computed from the evolutionary 
models of Bressan et al. ( 1993] ) for solar composition fol- 
lowing the interpolation algorithm described in Asiain et 
al. ( 1997 ). The algorithm considers, as input parameters, 
the T e ff and log g derived from the Stromgren photometric 
indices (Moon et al. 1985; Napiwotzki et al. 1993). 

The weakness in this procedure lies in the inability 
to take into account the effects of stellar rotation when 



deriving ages from photometry. Figueras & Blasi ( 1998 ), 
analysing a sample of main sequence B7-A4 stars, found 
that actual photometric ages increment by 30-50% in aver- 
age, if rotation is not considered. As individual corrections 
are not possible, and an important fraction of the 0-B9 
main-sequence stars are high rotators, this important sys- 
tematic trend has to be considered when deriving an age 
for the Gould Belt system. 

First, individual ages were derived for all the stars for 
which photometric data was available, without taking into 
account binarity, variability or photometric peculiarities. 
The age and relative error in age distributions computed 
for 2864 stars are presented in Fig. || and Fig. 0, respec- 
tively. With the aim of retaining as many as possible of the 
very young stars in our final sample, a careful treatment 
was followed to take into account the effects of binarity, 
duplicity or peculiarity in the age computation: 

— In the case of double or multiple stars for which 
only joint photometry is available (i.e., systems with 
p < 10" and AiJp < 3 m ), the computed age is greater 
than the real age (Trimble & Ostriker 1981). All dou- 



ble or multiple systems with computed ages smaller 
than 30 Myr have an actual age inside the interval t < 
30 Myr, so they were retained in the working sam- 
ple only if they have a reliable Hipparcos trigonomet- 
ric parallax [a^jir < 0.25), as the effects of duplicity 
render the photometric distance meaningless. Further- 
more, roughly 80% of these stars have t + a T < 60 
Myr, so there is a high probability that they belong to 
the Gould Belt. Double or multiple systems with an 
estimated age larger than 30 Myr were rejected. 
We rejected those stars with photometric or spectral 
peculiarities. The photometric indices accounting for 
temperature for Bp stars are bluer due to peculiari- 
ties, thus producing estimated ages smaller than real 
(Hauck |1975| ). Variable stars with AH p > 0.6 m were 
also rejected. 

As shown in Fig. ^, about 12% of the sample has a 
relative error in age larger than 100%. From them, 
those placed below the ZAMS (206 stars) are pecu- 
liar stars that are already rejected or very young stars. 
We checked that the stars above the ZAMS had a com- 
puted age smaller than about 30 Myr (116 stars). In 
any case, as these stars are expected to be very young 
and therefore of great importance to our study of the 
Gould Belt (they represent roughly 25% of the stars 
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Fig. 6. Age distribution of the sample of O and B stars 
(2864 stars). 
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Fig. 7. Distribution of the relative error in age of the O 
and B stars in the sample. 88% of the stars have a relative 
error in age smaller than 100%. 

Table 1. Averaged errors in the samples of O and B stars. 



Error 


Sample 


1 


Sample 2 


7j,r 


0.60 mas 




0.57 mas 


Jit) 


0.168 




0.163 


G fi a cos 5 


0.83 mas 


yr- 1 


0.81 mas yr -1 




0.70 mas 


yr- 1 


0.67 mas yr - 








3.44 km s" 1 



with r < 30 Myr), we decided to retain them in the 
final sample. 



2.4- Working samples 

Following the procedure described above, two samples of 
stars were formed: 
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Fig. 8. Distance distribution of the samples 1 (top) and 2 
(bottom) of O and B stars defined in the text. Filled his- 
tograms show the stars for which Hipparcos trigonometric 
distance was used (56.3% and 50.2% of the stars belonging 
to samples 1 and 2, respectively). 
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Fig. 9. Averaged errors in the three velocity components 
plotted against heliocentric distance. 

— Sample 1: Containing 3915 stars with known distance 
and proper motions. 

— Sample 2: A subsample of sample 1 containing 2272 
stars with known distance, radial velocity and proper 
motions. 

In Fig. U we show their distance distribution. Although 
the initial sample contains all the survey Hipparcos stars 
(complete up to V = 7.9), the lack of photometry and ra- 
dial velocity data reduces the completeness of the samples 
to about Vn m « 6.3, that is up to 150 pc for a B9V star, 
the faintest in our sample. 
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The averaged errors of several quantities are presented 
in Table |. In Fig. | we show the averaged errors in the 
three velocity components plotted against heliocentric dis- 
tance (computed taking into account the correlations be- 
tween the different variables provided by the Hipparcos 
Catalogue). 

The number of stars in samples 1 and 2 is reduced to 
2468 and 1789, respectively, when individual ages are re- 
quired. In Fig. [To] we show the position of these stars pro- 
jected on the X-Z plane (X positive towards the galactic 
center and Z towards the north galactic pole) , classified in 
three age groups (r < 30 Myr, 30 < r < 60 Myr and r > 
60 Myr). The Gould Belt is recognized as a tilted structure 
with regard to the galactic plane (Z — 0), mainly in the 
region with X < and Z < 0. This structure is clearly 
visible for stars younger than 60 Myr. The presence of 
some stars belonging to this structure in the interval 60 
< t < 90 Myr is fully justified by the large errors in the 
estimate of individual ages. 

3. The Gould Belt 

3.1. Structural parameters of the Gould Belt 

A classic problem in the study of the Gould Belt has been 
how to separate the stars belonging to this structure from 
t hose belonging to the g alactic belt. Stothers & Frogel 
( |1974|) and Taylor et al. fll987| ) proposed different algo- 
rithms based on the individual assignation of the stars 
to either belt. More recently, Cabrera-Caho et al. (1999; 
hereafter CEA) proposed least mean classification error 
decision criteria to separate those stars belonging to each 
belt. These criteria, based on the spatial distribution of 
the stars, run into an obstacle: the classification of stars 
lying in the overlap region between both belts. 

In this paper, we followed an alternative approach pro- 
posed by CTG. The method assumes that the belts form 
two great circles in the celestial sphere, with a star density 
which decreases with the angular distance from each equa- 
tor belt. The decrease in star density is assumed to follow 
a Gaussian law, the standard deviation being the angular 
halfwidth of the belt. Therefore, we suppose that the den- 
sity distribution of the sample in the celestial sphere can 
be written as: 



<j(l,b) = a G {l,b) + a s (l,b) 



(2) 



where o~g and cr g are the density distributions around the 
Gould Belt and the galactic belt equators, respectively. 

This decomposition allows us to derive several parame- 
ters for the Gould Belt: the spatial orientation with regard 
to the galactic plane (ic, ^g)> the fraction of stars belong- 
ing to the Gould Belt (q), and the angular halfwidth of 
each belt (£g ; £g)- Comeron (1992), using pre-Hipparcos 
data, showed that the q parameter describes in an ap- 
propriate way the main characteristics of the distribution 
of young stars, in particular the extent of the belts as a 
function of galactic longitude. 



The resolution procedure, based on the maximum like- 
lihood method, can be found in Comeron ( |1992| ). Addi- 
tionaly, an iterative procedure until convergence was im- 
plemented here to minimize the dependence of the final 
results on the departing values. As explained in Comeron 



(1992), the method requires the homogenous completeness 
of the sample throughout over the entire celestial sphere. 
Unfortunatelly, although Hipparcos is complete up to V = 
7.9, the absence of photometric measurements reduces this 
limit, so substantially reducing the distance limit of the 
intrinsically faint stars. As a good compromise between 
the constraints of completeness and the need for a statis- 
tically representative number of stars we considered only 
those stars that are brighter than V = 7.0 in sa mple 1. 
Numerical simulations, presented in Appendix A:, allowed 



us both to assess how these incompleteness effects could 
influence the results and to have an external evaluation of 
the errors on the derived structural parameters. 

The results obtained for the real sample are presented 
in Table § We see that the Gould Belt's structure is clearly 
detected in the subsamples of young stars with R < 600 
pc. From the top left panel in Fig. [To] (stars with r < 
30 Myr) we also see that this structure extends up to 
600 pc into the south galactic hemisphere and only up 
to 200-300 pc in the north. As young stars in the galactic 
plane reach distances greater than 1000 pc without having 
any substantial decrease in density, we can state that the 
distance cut off derived for this structure is real in this age 
interval and not a consequence of the incompleteness of 
our sample. An extent of about 600 pc is in agreement with 
Lindblad et al. (1997), who assumed that the prominent 



associations in the Gould Belt are within a distance of 700 
pc. 

The orientation parameters were found to be iq = 16- 
22° and = 275-295°. The orientation parameters are 
maintained up to the interval 60-90 Myr. For stars older 
than 90 Myr the method does not converge - large er- 
rors are found for ia and Oq, and q and halfwidths are 
undetermined -, so we conclude that the Gould Belt is 
no longer present. However, at this stage, it is mandatory 
to verify that the dissapearance of the structure for older 
stars is not a consequence of our observational constraints 
(in our sample, stars with r w 60 Myr have a limited dis- 
tance about 400 pc). Simulations show that distance cut 
off for these stars does not significantly disturb the de- 
termination of structural parameters. In other words, if a 
substantial number of stars with ages larger than 60 Myr 
would be present in the Gould Belt, our algorithm would 
be able to detect them from the number of stars available 
at present. The smaller ic value derived in the interval 
30-60 Myr compared with that obtained in the intervals 
r < 30 Myr and 60 < r < 90 Myr has no explanation 
at present (it cannot be attributed to the sample distance 
horizon). From the simulations, an uncertainty of 3-6° is 
expected for zq in these age intervals. In any case, our 
values for the orientation parameters are in good agree- 
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Fig. 10. Distribution on the X-Z galactic plane for the stars in the samples 1 (left) and 2 (right). At the top, the 
stars with an age less than 30 Myr; in the middle, those with an age between 30 Myr and 60 Myr; and at the bottom, 
those with an age larger than 60 Myr. 

Table 2. Structural parameters of the Gould Belt as a function of distance and age: inclination (ic), longitude of the 
ascending node (f^c), fraction of stars belonging to the Gould Belt (q), angular halfwidth of the Gould Belt (£q) an d 
the galactic belt (£ g ), and number of stars (N). Limiting visual apparent magnitude: 7.0. 



R (pc) 


iG n 


«G (°) q 


cg n 


£g (°) 


N 


t < 30 Myr 


R < 400 


2L2( 13 ) 


287.3(4.2) 0.60 


6.2 


22.5 


236 


R < 600 


19-9(1.6) 


282.8( 5 .2) 0.66 


7.0 


22.6 


300 


600 < R < 2000 


11-8(2.2) 


316.1(10.7) 






126 


30 < r < 60 Myr 


R < 400 


15-9(2.5) 


294.9( 6 .5) 0.64 


7.3 


23.7 


261 


R < 600 


15-5(2.6) 


293.4(6.5) 0.62 


7.2 


22.5 


297 


600 < 7? < 2000 


11-9(22.7) 


192.6(i64.i) 
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60 < r < 90 Myr 


R < 400 


22.3(2.i) 


276.1(4. 9 ) 0.44 
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25.8 
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22.1(4.1) 


276.7(4.1) 0.42 
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25.7 
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R < 400 


5-0(32.8) 


3 1 6.7(378.4) 
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3-4(33.8) 


319.3(569.2) 
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ment with those published in the literature. Lesh ( 1968| ), 
Stothers & Frogel ( |1974| ) and Westin ( |1985| ) reported val- 
ues of i G = 19-22° and fi G = 270-300°. CTG reported 
i G = 22.3° and fl G = 284. 5° from a sample of O-A0 stars. 
Recently, Guillout et al. ( |1998a[ ) analysed the sky distri- 
bution of X-ray emitting stars belonging to the RASS- 
Tycho sample and found a low galactic latitude feature - 
the Gould Belt - with an orientation with regard to the 
galactic plane of i G = 27.5 ± 1° and n G = 282 ± 3°. Even 
if we take the errors into account, this inclination value 
is not compatible with ours. A possible explanation for 
this discrepancy is the fact that Guillout et al.'s sample 
was restricted to stars that were not located at such great 
distances (R < 200 pc), where the most prominent struc- 
ture belonging to the Gould Belt - the Sco-Cen complex 
- defines a slightly higher slope than those located at a 
greater distance and in the opposite direction (mainly the 
Ori OBI association). CEA, after identifying those stars 
belonging to the Gould Belt, found i G = 17.5-18.3° and 
fie = 287-294°, which is in very close agreement with our 
results. 

For stars younger than 60 Myr, the fraction of stars 
belonging to the Gould Belt (q) were found to be 0.60-0.66. 
This value of q decreased to 0.42-0.44 when c onsidering 
st ars w ith an age between 60 and 90 Myr. From Appendi^ 



increase that can be evaluated in 30-40 Myr (Figueras & 
Blasi |1998|). 



A: we e xpect an uncertainty of « 0.15 in q. Therefore, its 



decrease for old stars must be real. CEA, classifying by 
spectral type groups, found q — 0.44 for 0-B2.5 stars with 
R < 1000 pc, and q = 0.36 for 0-B9.5 stars. 

The angular halfwidths of the belts were found to be 
£g = 6-8° and £ g = 22-26° (s ee Table p|) w ith an approxi- 
mated uncertainty of 5° (see Appendix A:). These values 
tend to increase smoothly when we go from young to not 
so young stars. This effect, if real, could be simply due 
to the geometrical effect produced by the decrease in the 
mean distance when older stars are considered. We found 
that the Gould Belt is narrower than the galactic belt, 
contrary to the result obtained by CTG, and in perfect 
agreement with the relationship 1:3 obtained by Stothers 
& Frogel ( 1974 ) studying the scale heights for both belts. 
The same trend is observed in the paper by CEA. 

3.2. The age of the Gould Belt 

On the basis of the spatial distribution analysed in the 
previous section, and using individual photometric ages, 
we estimate the belt to be younger than 60 Myr. As ex- 
plained before there are two major biases in the compu- 
tation of individual photometric ages that account for the 
presence of certain stars belonging to the Gould Belt in 
the age interval 60-90 Myr. First, due to the significant un- 
certainties in the age computation, we expect that some 
stars with actual ages smaller than 60 Myr are included 
in this interval. Second, and more importantly, the fact 
that we cannot include the effects of stellar rotation on 
the age computation scheme means there is a systematic 



The age estimates for the Gould Belt found in the lit- 



erature lie in the interval 20-90 Myr. Lesh (1968) found an 



age of 45 Myr when supposing two superposed stellar dis- 
tributions, and 90 Myr when considering only one homo- 
geneous population in expansion. Lindblad et al. ( 1973 ) 
suggested an age of 30-40 My r, the system being born 
in a spiral arm. Franco et al. ( |L988[ ) found an age of 60 
Myr for the Orion and Mo nocer os molecular cloud com- 
plexes. Tsioumis & Fricke ( |1979| ) and Comeron & Torra 
( 1991 ) reported 60 Myr and 70 Myr, respectively, from 
kinematic studies. From stellar individual age determina- 
tions using Stromgren photometry, Westin (1985) found 
an upper limit of 60 Myr. CTG divided their sample of 
O and B stars in spectral type subsamples, and reported 
a lower limit for age equal to the lifetime of a B4-type 
star (about 50 Myr). Guillout et al. (1998b), analysing 
the RASS-Tycho sample, reported the detection of a very 
young active late-type stellar population belonging to the 
Gould Belt. The X-ray luminosity distribution was com- 
patible with an age of 30-80 Myr for these stars. More 



recently, Moreno et al. (1999) derived an age of 20 Myr, 
assuming that the Gould Belt was formed by an expanding 
shell. 

A parallel determination of the Gou ld Belt's age from 



its kinematics will be given in Sect. 4.3 



4. Kinematics of young stars in the solar 
neighbourhood 

4-1- Fit of the kinematic model 

4.1.1. Equations for the linear model 

The Oort constants were derived using the first-order de- 
velopment of the systematic velocity field: 

V T = AR sin 21 cos 2 b 

+C R cos 21 cos 2 b + KR cos 2 b 

—Uq cos I cos b — Vq sin I cos b — W© sin b 



Rk hi cos b — AR cos 21 cos b + B R cos b 
—C R sin 21 cos b 
+Uq sin / — Vq cos I 



(3) 



(4) 



i?fc/ib = — AR sin 2Z sin 6 cos & 

— C R cos 21 sin b cos b — K R sin b cos b 
+Uq cos I sin b + Vq sin I sin b — Wq cos b 

(5) 

where I, b are the galactic coordinates, R the heliocentric 
distance in pc, V T the radial velocity in km s _1 and fj,\, /^b 
the proper motions in " yr _1 of each star. The constant k 
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= 4.741 km yr (s pc [/©, Vq and Wq are the com- 

ponents of the peculiar motion of the Sun in km s -1 with 
regard to the circular velocity and A, B, C and K are 
the Oort constants, linear combinations of the gradients 
of the systematic velocity. No systematic motions perpen- 
dicular to the galactic plane were considered other than 
that arising from the solar peculiar motion. As recently 
verified by Palous (1995), the E, D and H terms describ- 



ing this systematic motion do not improve the results and 
their values are of low significance. 

4.1.2. Resolution procedure 

A weighted least squares fit was performed to estimate 
the model parameters from Eqs. (||),((l|) and (||), taking 
the residual velocity of the star as a random error. Conse- 
quently, the weight of each equation was chosen as (Creze 
19731 ): 



Pi 



(6) 



, obs 



where cr bs are the individual observational errors in each 
velocity component of the star, calculated by taking into 
account the correlations between the different variables 
provided by Hipparcos Catalogue, and cr cos is the projec- 
tion of the cosmic velocity dispersion ellipsoid (<ttjj 0Vj 
aw) in the direction of the velocity component consid- 
ered. The detailed iterative procedure applied to simulta- 
neously derive the model parameters and the cosmic dis- 
persion for each of the subsamples considered is explained 
in Sect. I4.1.3L 



To check the quality of the least-squares fits we con- 
sidered the x 2 statistics for N — M degrees of freeedom, 
defined as: 



N 

E 

i=l 



[Vi - y(xi;ai, ...a M )Y 



(7) 



where the independent data (sky coordinates and 

distances), yi the dependent data (radial and tangential 
velocity components), N the number of equations and M 
number of parameters to be fitted. According to Press 
et al. ( 1992p , if the uncertainties (cosmic dispersion and 
observational errors) are well-estimated, the value of \ 2 
for a moderately good fit would be \ 2 rs N — M, with an 
uncertainty of y/2 (N - M). 

To eliminate the possible outliers present in the sample 
due to both the existence of high residual velocity stars 



(Royer 1997) or stars with unknown large observational 
errors, we rejected those equations with a residual larger 
than 3 times the root mean square residual of the fit (com- 
puted as y/\yi — y{xi; a±, ...clm)] 2 /N) and recomputed a 
new set of parameters. We checked than no more than 
10-15 stars from the total sample were rejected using this 
procedure, so ensuring the objectivity of the rejection cri- 
terion. 



Three distance intervals were considered: 100 < R < 
600 pc, 600 < R < 2000 pc and 100 < R < 2000 pc. The 
first two intervals allow us to determine the influence of 
the Gould Belt on the stellar kinematics, since as we have 
seen there is no evidence of the presence of this structure 
for R > 600 pc. The last interval provides a global vision of 
the kinematics of young stars in the solar neighbourhood. 
Stars with R < 100 pc were not considered because, on the 
one hand, they do not give us information about galactic 
rotation and, on the other hand, peculiar motions of these 
nearby stars (small errors, so large weight) can disturb the 
results. 

At this stage it is appropiate to clarify the physical 
meaning of the A, B, C and K terms in the two regions 
we want to analyse. In the region 600 < R < 2000 pc, not 
affected by the Gould Belt, these terms are expected to 
reflect the local shape of the rotation curve. Thus, assum- 
ing a smooth variation in space, they are properly called 
Oort constants and they account for the divergence (K), 
vorticity (B) and shear (A: azimuthally, C: radially) of 
the general velocity field of the galactic disk in the solar 
neighbourhood. Oiling & Merrifield ( 199§| ) distinguished 
between Oort constants for the local shape of the rota- 
tion curve and Oort functions when accounting for their 
variation with the galactic radius. On the contrary, in the 
region 100 < R < 600 pc, these first derivatives of the 
velocity field will include the peculiar velocity field asso- 
ciated with the Gould Belt. Thus, strictly speaking, they 
should not be called Oort constants and so we will refer 
to them as Oort parameters. 

Before analyzing the different fits performed there are 
two interrelated aspects that deserve special attention: the 
possible biases in the fitting parameters induced by the 
characteristics of the sample and our observational con- 
straints - irregular spatial distribution of the stars, in- 
completeness effects, biases in the availability of radial ve- 
locity data (already discused in Sect. 2. 2), etc. - and the 
systematic errors in the parameters induced by the pres- 
ence of observational errors in the right hand side of the 
Eqs. (||), ( 0) and (||), not considered in our least square 
fit. Creze ( 1970 ) derived approximate analytical correc- 
tions to evaluate the second aspect but, as he stated, we 
think that numerical experiments arc the best way to glob- 
ally evaluate and discuss all the effects present in our real 
sample. Furthermore, these simulations can give us new 
insights in a long standing problem: the discrepancies ap- 
pearing between the solutions obtained using radial ve- 
locity and proper motion equations. Following that, we 
performed several numerical simulations that are detailled 
The results are commented, in the next 



Appendix B 



paragraphs, along with the results from real data. 

In Table §, we compare the solutions obtained when 
considering only radial velocity data - Eq. (|^) -, proper 
motion data - Eqs. (Q) and (p|) solved simultaneously - or 
the combined solution - Eqs. (||), (Q) and (||) - for those 
stars in sample 2 with 600 < R < 2000 pc. As can be 
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Table 3. Oort constants and solar motion for stars in the sample 2 with 600 < R < 2000 pc. Units: A, B, C, K in 
km s _1 kpc -1 ; Uq, Vq, Wq, a in km s _1 . x 2 /-/V oq is the value of \ 2 divided by the number of equations. N is the 
number of stars. The errors in the fitted parameters were computed as of = cr 2 C~ , where C^ 1 = Covu/a 2 is the 
corresponding element in the covariance matrix and a the standard deviation of the measurements with unit weight 
(Linnik 1963). As expected, large errors in K and Wq are obtained when derived from proper motion data or radial 
velocity alone, respectively. A constant cosmic dispersion of (crjy, ay, &w) — (8,8,5) km s _1 was considered in all 
cases. 



Components 


A 


B 


C 


K 


Uq 


^0 


Wq 


a 




N 










600 < 7? < 2000 


pc 










V T 


H.7 ( i.d 




-0.7(i.i) 


—3.0(0.7) 


7.2(i.3) 


16.6(1.2) 


7-4(3.2) 


14.3 


2.64 


308 


Ml + Mb 


14.5 { o.9) 


-12.8(o.6) 


0-7(0.9) 


— 0.7(3.2) 


10.7(i.o) 


12-1(0.9) 


8-5(0.5) 


11.9 


2.05 


308 


V x + Ml + Mb 


13.0(0.7) 


— 12-7(0.8) 


0-0(0.7) 


—3.3(0.7) 


8-7(0.9) 


14.1(0.8) 


8-6(0.6) 


12.7 


2.39 


308 






600 < R < 2000 


pc excluding the reg 


ion with 200 </ < 250° 






Vx 


14.0(i. 6 ) 




-1-0(1.2) 


-1.4(i.o) 


10.1(1.7) 


18.5 ( i. 6 ) 


8.2(4.o) 


15.1 


2.91 


224 


Mi + Mb 


13.5(i.i) 


— 11-8(0.9) 


-1-1(1.4) 


-9.4(4.2) 


8.6(i.4) 


13.7(i.3) 


8-9(0.6) 


13.1 


1.19 


224 


V x + Ml + Mb 


13-5(1.0) 


-ll-4(i.o) 


— 1-1(0.9) 


—2-3(0.9) 


8-5(1.2) 


15.8(i.o) 


8-9(0.8) 


13.5 


1.25 


224 



seen, a difference about 2-3 km s -1 kpc -1 in the A value is 
present between the radial velocity and the proper motion 
solutions when solved separately. One effect that coul d 
reduce this discrepancy was pointed out by Creze ( 1970|) : 
a non-negligible standard error in distances produces an 
understimation of the A Oort constant when derived from 
radial velocities. Our numerical experiments (see Table 
Bl) indicate that this effect is less important than the 
bias induced by the cut in the observed distance (affected 
by errors), which produces a bias of about 1-1.5 km s _1 



crepancy on A vanishes (see Table |3|) . We also confirmed 
that when other special regions are removed, the discrep- 
ancy is maintained. Another important point we realize 
when studying these results is the variation in % 2 statis- 
tics when rejecting those stars in this region. For proper 
motion data, X 2 /A r C q decreases from 2.1 to 1.2. On the 
contrary, for radial velocity data there is a little increase 
from 2.6 to 2.9. For the combined resolution a decrease 
from 2.4 to 1.3 was found. It seems indicate that there is 
a better fit of the velocity field to the model when rejecting 



could ( ven be enlarged if the observational bias could be 
removed. 



kpc -1 in the opposite sense. From these simulations we this region. However, the value of x /N eq obtai ned is still 
conclu de that the difference present in the real sample lar ger than the derived in the simulations (see Appendix 

B:). This could be due by an underestimation in the error 
in the photometric distances and/or radial velocities for 
our stars. 



Recently, Feast et al. (1998), working with a sample 



of Hipparcos Cepheids, showed that the discrepancy ap- 
pearing between the A val ue derived from proper motions 
(Feast & Whitelock 1997 ) and that derived from radial 
velocities (Pont et al. 1994 ) disappears when the new Hip- 
parcos Cepheid distance scale is considered in the last 
equations. In our case, we verified that an overestimation 
in our photometric distances by a factor of 20% - e.g. as- 
suming a stellar rotation effect in the absolute mag nitude 
derivation (Lamers et al. 1997 , Domingo & Figueras 1999 ) 
- can account only for a difference of 1-2 km s _1 kpc -1 
between both solutions. 



Recently, Palous (1998) used Hipparcos data to esti- 



Lindblad et al. (1997) attributed this discrepancy to 
the irregular distribution of stars and stellar groups and 
also to a possible non-linearity in the velocity field. As ex- 



plained in Appendix B: our simulations take into account 



the irregular spatial distribution of real stars and we saw 
there that the discrepancy did not appear. More promiss- 
ing is the hypothesis that the discrepancies come from 
departures of some stellar groups from the adopted linear 
model. In our sample, we verified that when the stars in 
the region 200 < I < 250° (partly composed by stars of the 
Ori OBI and Col 121 associations) are removed, the dis- 



mate the second-order terms in the expansion of the ve- 
locity field of young stars around the Sun. He concluded 
that these terms are always of low significance and do not 
significantly alter the values obtained for the first order 
derivatives (Oort constants). It cannot be ruled out that 
higher order terms might account for the irregularities 
observed, particularly the discrepancies appearing when 
some regions in I are omitted. However, the small num- 
ber of stars available at large distances does not allow us 
to avoid the large correlations between variables when sec- 
ond order terms are included. We believe that the peculiar 
motions in certain specific regions are largely responsible 
for the discrepancies observed in the derivation of the A 
Oort constant. 

Likewise, we confirmed that the discrepancy on Vq 
does not disappear when eliminating the stars of the Orion 
region, neither when eliminating stars from other regions. 
Again, from our simulations we discarded the irregular dis- 
tribution of the stars in the X-Y plane as an explanation 
for this discrepancy, so again the most probably explana- 
tion is the departure from the adopted model. Finaly we 
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would to comment that the difference present in the de- 
rived K value is a simple consequence of the fact that this 
parameter is poorly determined when only proper motion 
equations of stars with small galactic latitude are consid- 
ered. 

Following we present the correlation matrices obtained 
for the radial velocities, proper motions and combined 
solutions given in Table |: 



Table 4. Standard deviation of the observational errors 
and cosmic dispersion for several subsamplcs divided in 
age intervals in the distance interval 100 < R < 600 pc, 
both expressed in the galactic heliocentric coordinates sys- 
tem. er cos = (afj +<jy + cr^) 1 / 2 : (1): this work, (2): Wielen 
(1977) analytical approximation er(t)" = er™ + C„t with 



n = 2 and C v = 6 • 10~ 7 (km s" 1 ) 2 
s _1 . Units: km s _1 . 



yr 



_1 and (To = 10 km 



Radial velocity solution. 










Age 
(Myr) 




(ft/, °V, &w)cos 


^"cos 

(1) 


^COS 

(2) 


A C K 


Uq 


Vq 


W @ 




- 30 


(3.5, 3.0, 2.3) 


(7.9, 7.2, 4.3) 


11.5 


10.4 


1.00 -0.12 -0.07 - 


-0.18 


-0.15 


-0.12 


A 


30 - 60 


(3.1, 3.4, 1.8) 


(6.2, 7.5, 4.4) 


10.7 


11.3 


1.00 -0.03 


0.05 


-0.04 


-0.00 


C 


60 - 90 


(3.2, 3.2, 2.0) 


(7.5, 8.8, 4.5) 


12.4 


12.0 


1.00 - 


-0.11 


-0.06 


-0.09 


K 


90 - 120 


(3.2, 3.1, 2.1) 


(10.9, 9.6, 6.4) 


15.9 


12.8 




1.00 


-0.08 


-0.04 


Uq 


> 120 


(3.2, 3.2, 2.9) 


(10.8,10.0, 5.5) 


15.7 


13.8 






1.00 


-0.05 
1.00 


Vq 
Wq 













ii). Proper motion solution: 



A B C 


K 


u & 


v & 


Wq 




1.00 -0.14 0.11 - 


-0.06 


0.05 


0.07 


0.06 


A 


1.00 0.14 


0.03 


-0.02 


-0.15 


-0.01 


B 


1.00 


0.04 


-0.22 


0.10 


0.01 


C 




1.00 


-0.04 


-0.08 


0.24 


K 






1.00 


0.10 


0.04 


Uq 








1.00 


0.05 


Vq 










1.00 


Wq 


in). Combined solution. 












A B C 


K 


Uq 


Vq 


Wq 




1.00 -0.09 0.02 - 


-0.07 


0.07 


0.08 


0.01 


A 


1.00 0.11 


0.01 


-0.02 


-0.09 


-0.01 


B 


1.00 


0.02 


-0.10 


0.02 


-0.01 


C 




1.00 


-0.08 


-0.07 


0.01 


K 






1.00 


0.01 


0.01 


Uq 








1.00 


0.05 


Vq 










1.00 


Wq 



the residual analysis of the combined solution. Henceforth, 
we will work hereafter with the combined solution solving 
Eqs. (||), (Q) and (||) simultaneously so taking into account 
all the information available. 



4.1.3. Cosmic dispersion 

Our residuals in Eqs. (||), ([|) and (||) include the observa- 
tional errors, the residual velocity of the stars (cosmic dis- 
persion) and possible departures from the adopted linear 
model. An optimal solution to characterize and separate 
the contribution of the cosmic dispersion may be to adopt 



As can be seen, the correlations are small in all cases, 
constating that they are not responsible of the discrepan- 
cies between proper motions and radial velocities. We also 
realized that the combined solution presents the smallest 
correlations. 

The information that can be derived from the \ 2 statis- 
tics provides us further arguments to favour the combined 
solution in the analysis performed in the following sec- 
tions. As can be seen in Table ||, the worst values for the 
fraction x 2 /N cq are obtained in the radial velocity solu- 
tion, effect that is not reflected in our simulations. Two 
aspects can contribute to this fact: an underestimation of 
the observational errors in radial velocities - which can 
not be rule out since this parameter is very difficult to ob- 
tain for very hot stars - or an error in the adoption of the 
shape and size of the velocity dispersion ellipsoid. Both as- 
pects will be discussed hereafter but the only way to fully 
analyze the cosmic dispersion velocity ellipsoid is through 



it from independent works such as Wielen (1977), Lacey 
( |1991D , Asiain et al. ( |l999h| ), among others, who analit- 
ically or empirically evaluated the increase of the stellar 
velocity dispersion with age (disk heating). Discrepancies 
are present in the first 100-150 Myr (our working domain), 
where the population is not relaxed. 

Working in the opposite sense, we can assume that the 
linear model adopted is approximately correct and that 
the observational errors are well estimated. With this hy- 
pothesis, and considering the interval 100 < R < 600 pc 
to minimize both the possible departures from the lineal 
model and the observational errors, we solved the com- 
bined solution for different age intervals using an iterative 
process up to convergence. In all cases one interation is 
sufficient and the results do not depend on the adopted 
initial values. The results are presented in Table |J 

We must point out that the values obtained for the 
Oort parameters and the solar motion components are 
practically independent of the choice of the cosmic dis- 
persion values, the differences in these parameters being 
always smaller than 0.5 km s _1 kpc -1 or 0.5 km s _1 re- 
spectively. Conversely, the adopted values for the cosmic 
dispersion will directly modify the \ 2 statistic. With this 
in mind, and the fact that the values derived for the cos- 
mic dispersion are coherent with Wielen's work (see Ta- 
ble U), we proceded to solve combined solutions at dif- 
ferent age intervals using the above explained iterative 
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process for the interval 100 < R < 600 pc. For the inter- 
val 600 < R < 2000 pc and 100 < R < 2000 pc we use 
as cosmic dispersion the values obtained for 100 < R < 
600 pc, that is, assuming isothermality. This process give 
us the oportunity to use the \ 2 statistic to evaluate the 
over/underestimation of the observational errors or the 
departure from the adopted linear model. The results arc 
summarized in Table || and discussed in the following sec- 
tions. 

4-2. Large scale outline of the local galactic kinematics 

An initial overview of the kinematics of young stars in the 
solar vicinity is given when considering all the stars in the 
interval 100 < R < 2000 pc. The solar motion relative to 
the stellar group considered was found to be: 

(U Q ,V Q ,W Q ) = (11.0, 12.9, 6.8) ±(0.2, 0.2, 0.1) km s _1 (8) 

On the other hand, we found the following values of the 
Oort constants, which are dominated by a pure differential 
galactic rotation: 

A = 11.8 ±0.4 km s" 1 kpc" 1 

B = -12.3 ±0.4 kms^kpc" 1 

C = 0.4 ±0.4 kms" 1 kpc" 1 

K = -2.0 ±0.4 kms^kpc -1 (9) 

The kinematic distortion produced by the Gould Belt 
in the solar neighbourhood can be removed when consid- 
ering only those stars with 600 < R < 2000 pc. Then, we 
found a solar motion: 

(U Q ,V G ,W Q ) = (9.0, 13.4, 8.3) ±(0.8, 0.7, 0.5) kms _1 (10) 

in perfect agreement with the classic value, the changes 
with respect to the former solution being partly produced 
by the p resence of mov ing groups am ong young stars (Asi- 
ain et al. 1999a ). From Appendix B: we expect these values 
to be slightly underestimated due to a bias of « —(0.3-0.4) 
km s" 1 . The Oort constants were found to be: 



.4 
B 
C 
K 



13.0 ±0.7 kms" 1 kpc" 1 
-12.1 ±0.7 kms -1 kpc" 1 
0.5 ±0.8 kms" 1 kpc" 1 
-2.9 ±0.6 kms -1 kpc" 1 



(11) 



From Appendix B:, a significant bias is only expected for 



B, with an underestimation (in absolute value) of about 
0.8 km s _1 kpc -1 . Therefore, its value might be B « 
-12.9 ±0.7 km s" 1 kpc -1 . These A and B Oort constant 
values are in good agreement with the results obtained by 



Lindblad et al. (1997). Using Hipparcos data they found 
A = 13.7 ± 1.0 km s" 1 kpc" 1 and B = -13.6 ± 0.8 km 
s" 1 kpc" 1 from a sample of O and B stars with R < 2000 



km s _1 kpc" 1 from a sample of Cepheid stars (at dis- 
tances up to 5000 pc) with Hipparcos proper motions and 
distance calibration. Using a similar sample (also with an 
Hipparcos distance calibration), Feast et al. ( [1998 ) found 
A = 15.1 ± 0.3 km s" 1 kpc" 1 from radial velocities. This 
tendency to obtain lower values for the A Oort constant 
when the distance horizon of the sample is approached is 
confirmed in the results presented in Table [|. Using stars 
with 100 < R < 600 pc not belonging to the Gould Belt 
(age larger than 90 Myr), we found an A constant of 11.9 
± 2.0 km s" 1 kpc" 1 , roughly 1 km s" 1 kpc -1 less than 
that obtained for all the stars with 600 < R < 2000 pc. 
When comparing this result with the A = 11.3 ±1.1 km 



s 1 kpc 1 derived by Hanson (1987) using proper motions 



for approximately 60 000 nearby faint stars (R < 1000 pc 
and a photographic magnitude 16 < m pg < 17) from the 
Lick Northern Proper Motion (NPM) program, we found 
a good coherence, although his sample was compose d bas i- 
cally of F2-K0 stars. Recently, Oiling & Merrifield ( [1998D , 
using a mass model which includes the interstellar gas 
component, derived the variation in the Oort constants 
as a function of the galactocentric distance. The authors 
explained the discrepancies between the A values derived 



by Hanson (1987) and Feast et al. (1998) as being pro- 
duced by the different mean galactocentric distances of 
the two samples. However, the authors admitted a poten- 
tial source of error in their analysis due to the assump- 
tion of azimuthal symmetry in the orbital structure of the 
Galaxy. 

Our nearly null values of C and K constants when con- 
sidering all the stars in the distance intervals 100 < R < 
2000 pc and 600 < R < 2000 pc are in good agreement 
with a pure differential galactic rotation. CTG found a 
value of K w — (1-2) km s _1 kpc -1 for stars with R < 
1500 pc, very similar to ours. Nevertheless, for the C con- 
stant they found a clearly negative value for B6-A0 stars: 
C = —8.8 ± 1.1 km s" 1 kpc -1 . This negative value of C 
seems to be corroborated by Mestres ( 1996] ), who found 
1 kpc -1 for those stars with 400 



C = -4.8 ± 1.2 km s 

< R < 1500 pc. More recently, Lindblad et al. (|L997j) 

- 1 kpc" 1 and K = -1.1 ± 0.8 



found C = 0.8 ± 1.1 km s" 



pc outside the Gould Belt. Feast & Whitelock ( |1997| ) re- 
ported A = 14.8 ±0.8 km s" 1 kpc" 1 and B = -12.4 ±0.6 



km s kpc" from their sample of O and B Hipparcos 
stars outside the Gould Belt. This last result is in very 
good agreement with ours. 

To evaluate the goodness of our fit we can look at the 
X 2 statistics. We found values 1.9-2.3 for \ 2 /N cq , depend- 
ing on the age interval considered. For a moderately good 
fit we would expect a value x 2 /N c(l 1. We think the 
difference to be produced by an underestimation in the 
errors in the photometric distances and radial velocities 
for far stars. 

This overall view is changed considerably when we di- 
vide our sample of stars into age and distance groups. The 
most characteristic change when we study the system of 
the nearest and youngest stars is the kinematic signature 
of the Gould Belt, i.e. the appearence of a non-null value 
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Table 5. Oort constants and residual solar motion as a function of distance and age. Units: Age in Myr; A, B, C, K 
in km s _1 kpc -1 ; Uq, V©, Wq, a in km s -1 . x 2 /N cq is the value of x 2 divided by the number of equations. N is the 
number of stars (sample 1 + sample 2). 





A 




C 


K 


u & 








A 1 1 c q 


N 










100 < R < 600 


pc 










0-30 


5.7(i.4) 


-20.7(1.4) 


5.2(i. 4) 


7.1(1.4) 


8-1(0.5) 


14-5(0.4) 


6-4(0.3) 


6.02 


1.06 


361 + 289 


- fin 

OU \J\J 


7-6(1.5) 




y -°(1.6) 


4.0 ( i. 7) 


11-6(0.4) 




' -^(0.3) 








< 60 


6.3(i i) 


-18.5(i.o) 


5.9(i.i) 


5.1(i.i) 


9-8(0.3) 


14 4m i\ 


6 9m o\ 


6.12 


1.01 


720 + 555 


60 - 90 


10.5(2.i) 


— 13.6(2 o) 


5.9(2.i) 


-5.4(2.3) 


12.4(„.5) 


lo.O/n ^^^ 

v-^O.Oj 


O.of n /i \ 


6.72 


1.14 


245 + 183 


> 60 


11.8 n « 


-11.0(1.4) 


-0.9(i. S ) 


0.0(1.7) 


ion 
J-^.U(Q.4) 


t O . Z l(\ A\ 
" ■ 1 U . 4 J 




8.07 


1.09 


932 + 654 


> 90 


11 9t9 nl 


-9 4m «i 


— 4 Qio ri\ 


1.9(2.2) 


11-8(0.4) 


1 Z 9 m .1 1 


6 6m ^> 


8.30 


1.09 


687 + 471 


All 


8. Off) R) 


-14.2(o.7) 


1.5(0.8) 


0.5(0.9) 


11-2(0.2) 


13.0(o.2) 


6.7(o.i) 


7.24 


1.06 


2970 + 1596 










600 


< R < 2000 


pc 










- 30 


13.3(o.7) 


-11-7(0.7) 


-0.3(0.7) 


—2.6(0.7) 


8-0(0.8) 


12.9(0.8) 


7.9 { o.5) 


10.39 


1.95 


285 + 204 


30 - 60 


9.1(1.7) 


-10.8(i. 8 ) 


-3.4(i. 8) 


-0-8(1.8) 


15-0(1.8) 


10.4(i.7) 


8.9(i.i) 


11.57 


2.29 


81 + 56 


< 60 


12-7(0.6) 


-11-7(0.7) 


-0.6(o.7) 


— 2.5(0.6) 


9-1(0.8) 


12.5(o.7) 


8-0(0.5) 


10.91 


2.05 


366 + 260 


All 


13.0(0.7) 


-12.1(0.7) 


0-5(0.8) 


—2-9(0.6) 


9.0(0.8) 


13.4(„.7) 


8-3(0.5) 


11.77 


1.87 


449 + 308 










100 


< R < 2000 


pc 










- 30 


12.9(o.6) 


-13.0(0.6) 


0.5(0.6) 


— 1-7(0.5) 


8-6(0.4) 


13.4(„.4) 


6-7(0.3) 


7.85 


1.51 


646 + 493 


30 - 60 


9-6(i.o) 


— 13-2(1.0) 


2.0 (1 . ) 


-0.2 ( i.o) 


11-7(0.4) 


13.4(0.5) 


7-4(o.3) 


6.85 


1.19 


440 + 322 


< 60 


12.0(0.5) 


-13.0(0.5) 


0.7(0.5) 


-1-5(0.5) 


10-0(0.3) 


13.4(o.3) 


7-1(0.2) 


7.52 


1.34 


1086 + 815 


> 60 


11.1(1.4) 


-12.2(i.3) 


— 1-8(1.3) 


-5-1(1.5) 


11-8(0.4) 


13.2(o.4) 


6.8(o.2) 


8.45 


1.19 


981 + 676 


All 


11-8(0.4) 


— 12.3(0.4) 


0.4(0.4) 


— 2.0(0.4) 


11-0(0.2) 


12.9(0.2) 


6.8(o.i) 


7.77 


1.17 


3419 + 1904 



of the K Oort constant and the peculiar behaviour of the 
other Oort constants. 



4-3. Local irregularities: the kinematic effects of the 
Gould Belt 

To study the kinematic characteristics of the Gould Belt, 
we present in Table || the results for those stars in our 
sample with 100 < R < 600 pc. In Fig. |ll|, we show the 
variation of the Oort parameters with age. In general, we 
observe a marked increase in A and B values with age, and 
a decrease in C and K, according to the results obtained 



by Torra et al. (^997J). 

A non-pure differential galactic rotation was found for 
the youngest group of stars, with A — 5.7 ± 1.4 km s _1 
kpc^ 1 , B = — 20.7 ± 1.4 km s _1 kpc -1 and non-null values 
for C and K. The tendency to obtain small A and B values 
for the youngest group is in perfect agreement with the 
results obtained by Lindblad et al. (1997) and Torra et 
al. (1997), although these autors even report negative A 



values in the combined solution. In agreement with these 
two studies, we confirmed that significant differences do 
appear in the solutions using only radial velocities or only 
proper motion equations (A — 0.7 ± 2.9 km s _1 kpc -1 
and A — 6.4 ± 1.4 km s" 1 kpc -1 , respectively), and as 
discussed in the previous section, a possible explanation 
could be the departure of some stellar groups from the 
adopted linear model. When we considered not so young 
stars, A and B are closer to classic values (A m 12 km s -1 
kpc -1 and B pa —9 km s _1 kpc -1 ). From Appendix B: 



produced by any systematic bias. In this distance interval 
we only expect an underestimation in A of about 0.5 km 
s" 1 kpc -1 . 

For those stars with t < 60 Myr we found a clear 
positive value of K-term: K = 7.1 ± 1.4 km s -1 kpc -1 
for t < 30 Myr, and K = 4.0 ± 1.7 km s -1 kpc -1 for 30 

< t < 60 Myr. On the other hand, for stars older than 
60 Myr slightly negative values of K were found - nearly 
compatible with a null value. Positive values of K were 
not found when we considered the distance interval 600 

< R < 2000 pc, independently of the age interval. 

In the interval 100 < R < 600 pc the obtained values 
for Y 2 /N Pn were similar to those derived from the simu- 



we can see that these differences are not expected to be 



lations (see Appendix B: ), around 1.0. From the coherent 
values of the cosmic dispersion obtained, we can conclude 
that the velocity field of our stars fits the lineal model 
proposed and the errors are well estimated. 

The variations in Oort parameters as a function of 
age allow us to infer an estimation of the age of the 
Gould Belt. As we have seen, when considering stars with 
100 < R < 600 pc a nearly pure differential galactic rota- 
tion was only found for stars with an age greater than 90 
Myr. In the age interval 60-90 Myr, rather low values of 
A and B and a high value of C were still derived. So, we 
conclude that the age of the Gould Belt derived from the 
kinematic behaviour of the stars is in perfect agreement 
with that derived in Sect. |3.2| from the analysis of their 
spatial distribution. 

In Fig. [l2] we show the residual space velocity vectors 
for each star projected on the galactic plane, after sub- 
tracting the solar motion and the galactic rotation found 
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Fig. 12. Residual heliocentric space velocity vectors projected on the galactic plane for O and B stars with an age 
r < 30 Myr (top left), 30 < t < 60 Myr (top right), 60 < t < 90 Myr (bottom left) and 90 < t < 120 Myr (bottom 
right). 



Table 6. Oort parameters and residual solar motion for stars with 100 < R < 600 pc and r < 30 Myr when excluding 
those stars belonging to the complexes Sco-Cen and Ori OBI, following the member lists provided by Brown et al. 
( |1994| ; Ori OBI association) and de Zeeuw et al. (1999; Sco-Cen complex). Units: A, B, C, K in km s _1 kpc -1 ; U ( 



0) 



Vq, Wq, a in km s 1 . x 2 /-/V eq is the value of \ 2 divided by the number of equations. N is the number of stars (sample 
1 + sample 2). 



Excluded 


A 


B 


C 


K 


U e 


Vq 




a 




JV 


None 


5.7(i.4) 


-20.7(i.4) 


5.2(i.4) 


7.1(1.4) 


8-1(0.5) 


14.5(„.4) 


6-4(0.3) 


6.02 


1.06 


361 + 289 


Sco-Cen 


6.9(i. 6) 


-19.7(i.6) 


4-7(1.6) 


5.8(i.6) 


8-5(0.6) 


13.9(0.5) 


6-2(0.3) 


6.39 


1.18 


305 + 238 


Ori OBI 


6-1(1.6) 


-20.7(i.6) 


5.3(i.6) 


7-3(1.6) 


8-0(0.5) 


14.6(0.4) 


6-6(0.3) 


6.14 


1.10 


315 + 251 


Both complexes 


7-2(1.8) 


-19.7(i.8) 


4.9 ( i.9) 


6.0(i.9) 


8-4(0.0) 


14.0(o.6) 


6-4(0.3) 


6.60 


1.26 


258 + 200 
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Fig. 13. Heliocentric space velocity vectors (a,c) and residual velocity vectors (b,d) projected on the galactic plane 
(X-Y) and the meridional plane (X-Z) for stars younger than 30 Myr belonging to the complexes Sco-Cen and Ori 
OBI. 



o 




40 60 
Age (Myr) 



100 



Fig. 11. Variation of Oort parameters plotted against the 
age for stars within 100 < R < 600 pc. Units are in km 
s _1 kpc" 1 . 

in the resolution for all the stars in the distance interval 
600 < B < 2000 pc, classified in different age groups. From 



this figure, it is evident that to model the expansion of the 
Goul d Belt system as an expansion from a point (Olano 
[L982D or a line oriented in the direction I = 45° — > 225° 
(CTG) is not a good approximation. In the same figure we 
have identified certain clumps associated with Lac OBI, 
Cep OB2, Cas-Tau, Per OB2, Col 121, Vel OB2, Tr 10 and 
Sco-Cen. All these associations, exce pt Cep OB2, p ossibl y 
belong to the Gould Belt (Comeron |b99i| Poppel |1997|) . 
In the youngest group (top left panel: stars with r < 30 
Myr) we observe the clear residual motion of the Sco-Cen 
complex. This complex is composed of three associations 
- Upper Scorpius, US; Upper Centaurus Lupus, UCL; and 
Lower Centaurus Crux, LCC - at a distance of about 120- 
145 pc, in the region of positive X and negative Y. The 
Ori OBI association shows a smaller mean residual mo- 
tion, as will be discussed later. It is composed of several 
subgroups, situated in a range of distances of about 340- 
510 pc, in the direction of I w 200-210°. 

Sco-Cen and Ori OBI are the two main complexes 
in the Gould Belt, and it is particularly interesting to 
study their motion and influence on the velocity field sep- 
arately. This study might clarify whether the Gould Belt 
is a casual arrangement of OB associations or a struc- 
ture with a common origin. In Fig. [l^ we show the he- 
liocentric velocity field and the residual velocity field - 
computed as in Fig. [T^ of the stars belonging to the 
Sco-Cen and Ori OBI complexes on the X-Y and X-Z 
galactic planes. To select these stars, we used the member 



lists provided by Brown et al. (1994; Ori OBI association) 
and de Zeeuw et al. (1999; Sco-Cen complex). It should 
be noted that we have detected - especially in the case of 
Sco-Cen - some additional stars whose location and kine- 
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matics are compatible with membership to these associa- 
tions. We confirmed a high residual velocity field for Sco- 
Cen, which is moving away from the Sun: (U, V, W) les = 
(4.1, -6.7, 2.2) km s" 1 , (v T ,v h v h ) re J£= (7.1,-3.2,0.6) km 
s -1 . In the case of the Ori OBI association, the mean 
residual motion found was smaller and practically null in 
the radial direction: (U, V, W) ICS — (—2.7, 3.1, 3.6) km s _1 , 
(v T , v\, «b)rcs = (0.2, —3.9,3.8) km s" 1 . These effects were 
confirmed when deriving the Oort parameters for those 
stars with 100 < R < 600 pc and younger than 30 Myr 
not identified as members of these complexes (see Table 
^). When Ori OBI was excluded from the calculation of 
Oort parameters, the value of the K-tevm was found to 
rise only by 0.2 km s _1 kpc -1 . In contrast, when Sco-Cen 
was excluded a decrease of 1.3 km s _1 kpc -1 in the K 
value was found. In both cases the A, B and C Oort pa- 
rameters change less than about 1 km s -1 kpc" 1 . When 
both complexes were eliminated the parameters obtained 
were very similar to those obtained earlier (all the changes 
are within the error bars). 

We conclude that these associations are not the only 
responsible for the peculiar kinematics observed for the 
youngest stars in the solar neighbourhood, attributed to 
the Gould Belt. So, other nearby associations and field 
stars belonging to the Gould Belt have a great significance 
in the determination of the Oort parameters. 

An attempt was made at analysing the expansion of 
the system as a function of distance. As a first step, in Fig. 
|l4| we present the variation of the KR product as a func- 
tion of heliocentric distance for stars with r < 60 Myr. The 
radial expansion disminishes rapidly with increasing dis- 
tance (for R < 250 pc) and it does not extend further than 
400 pc. At distances larger than 300 pc, only Per OB2 has 
a mean residual motion away from the Sun. As discussed 
above, Ori OBI has an almost null radial residual motion. 
Even consi derin g the solar motion proposed by Dehnen 
& Binney fll99fl) - (U Q ,V Q ,W Q ) = (10.00,5.25,7.17) km 
s" 1 - and the galac tic rotation curve obtained by Feast 
& Whitelock ( 1997 ), we obtain a small residual motion 
for this agreggate of (U,V,W) ies = (-1.2,-2.8,2.1) km 



10 



(« r ,VL,Vb)r 



(1.7,1.9,2.8) km s" 1 . To analyse in 



more detail the expansion in the R < 300 pc region, we 
show in Fig. [l5] the residual velocity vectors projected on 
the galactic plane for stars with r < 60 Myr (for clarity 
we divided the sample in two different coronae: R < 150 



pc and 150 < R < 300 pc) together with Olano's ( [1982[ ) 
Lindblad Ring - with a center placed at R = 166 pc from 
the Sun in the direction I = 131°, and semiaxes of 364 pc 
and 211 pc - and the position of the center of the Gould 
Belt proposed by Comeron & Torra (1991) - R = 80 pc, 
I = 146°. In the first quadrant there is a lack of stars, 
partially produced by the near high extinction structures 
related to the Ophiuchus-Aquila complex (Vergely et al. 



1997), which does not allow us to study the residual ve- 



er 




100 200 300 400 
Heliocentric distance (pc) 



500 



600 



Fig. 14. Variation of KR against the heliocentric distance 
for stars with an age smaller than 60 Myr. 



locity field in this region. The Cas-Tau complex is moving 
in the direction I ps 240°. So, as shown in Figs. [l2| and 
p"5| and by Lindblad et al. (1997), the third and fourth 



galactic quadrants contain the most significant structures 
accounting for the expansion of the Gould Belt. 

Whereas the fourth quadrant contains the extensively 
studied Sco-Cen agreggate, a well-defined concentration 
of O and B stars with two different residual motions is 
present in the region 225 I 285°, mainly with dis- 
tances in the interval 100 Ss R Si 300 pc and ages between 
30 and 60 Myr. To analyse these streams in detail we 
present in Fig. [l6| the distribution of these stars i n the 
U-V plane, where a kernel estimator (Silverman 1986 ) 
was used to indicate the isocontours. Over this figure, 
the mean (U, V) heliocentric velocity components of the 
open clusters present in this region derived from Hippar- 
cos data by Robichon et al. ( 1999] ) and the new kine matic 
structures recently identified by Platais et al. ( |1998|) were 
also plotted. The stream placed at (U,V) = (—12,-23) 
km s" 1 shares the motion of the a Car (= HIP 45080) 
cluster and IC 2602 open cluster, and it may also be re- 
lated to the Pleiades moving group substructures found 
by Asiain et al. ( |l999a ). The concentration observed at 
(U, V) = (-28, -20) km s" 1 is associated with NGC 2451 
A and Tr 10. The nature of NGC 2451 has been under 
discussion for some time. According to Roser & Bastian 
(|1994D, NGC 2451 can be divided into two different enti- 



= 4.741-Rjti] cos b and v b = 4.741i?/i b 



ties. These authors named the closest of these the Puppis 
Moving Group (PMG), the center of its distribution be- 
ing clearly offset from the nucleus of NGC 2451 by 1°. 
The dista nce to PMG was found to be 220 pc. Carrier et 
al. ( 1999] ) also found two entities, at 198 and 358 pc re- 
spectively, from Geneva photometry and Hipparcos data. 
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Fig. 15. Residual heliocentric space velocity vectors projected on the galactic plane for O and B stars with an age 
r < 60 Myr and < R < 150 pc (left) and 150 < R < 300 pc (right). It is also drawn the Olano's ( |1982j) Lind blad 
Ring (with a cross in its center) and the center of the Gould Belt (square) obtained by Comeron & Torra (1991). GC 
and GR indicate the galactic center and the galactic rotation directions respectively. 



Two well-defined peaks in parallax seem to reinforce the 
nature of open cluster for these two entities, although the 
most distant is difficult to distinguish from the field stars 
because both the parallax and proper motion of its stars 
are close to those of field stars. On the other hand, Tr 10 
was identified as an intermediate age OB association by 
de Zeeuw et al. (199£), who found 23 members spread over 
~ 8° in the sky. Robichon et al. ( |1999| ) found 9 Hippar- 
cos members of the cluster. The distance derived in both 
papers is the same (365 pc). The U velocity component 
found for the association is the same as for the cluster 
(U = —27.3 km s _1 ), but there is a difference of 4 km s _1 
in the V component (Tas SO c = —17.8 km s _1 , V c \ u = —21.9 
km s _1 ). 

We found that only 7 stars in this region (225 < I < 
285°, 100 < R < 300 pc and 30 < t < 60 Myr) have been 
identified by the above mentioned authors as members 
of open clust ers or associations (including the list of de 
Zeeuw et al. 



1999) 



so a large number of the remaining 
stars share the motion of these clusters and associations 
and are spread over a large region, as observed in Fig. |l5| . 
Finally, although IC 2391 and HR 3661 (= HIP 45189) 
seem to be isolated in Fig. [l6], stars sharing the motion 
of these clusters were detected in our sample when the 
age interval was changed to 60 < t < 90 Myr. This is in 
agreement with the estimated age for HR 3661 - 100 Myr 
(Platais et al. |1998[) - but n ot with the age of IC 2391 - 



30 Myr (Stauffer et al. |i997| ) . 

Further work will be necessary to confirm the existence 
of these streams and to deal with their origin in the con- 
text of the various models proposed for the Gould Belt. As 



Table 7. Galactic coordinates and heliocentric velocity 
components of the clusters shown in Fig. [l6[ Units: I, b in 
degrees; R in pc; U, V and W in km s _1 . 



Cluster 


I 


b 


R 


U 


V 


W 


1. IC 2602^ 


289.6 


-4.9 


152 


-8 


-20 


-0 


2. a Car 1 


277.7 


-7.6 


132 


-11 


-24 


-4 


3. NGC 2232 2 


214.3 


-7.7 


325 


-16 


-12 


-11 


4. NGC 2516 2 


273.9 


-15.9 


346 


-17 


-24 


-4 


5. HR 3661 1 


266.9 


3.4 


174 


-22 


-15 


-6 


6. IC 2391 2 


270.4 


-6.9 


146 


-23 


-14 


-7 


7. Tr 10 2 


262.8 


0.6 


365 


-27 


-22 


-10 


8. NGC 2451 2 


252.4 


-6.8 


189 


-29 


-20 


-14 



1 Platais et al. (|1998|) ~ 

2 Robichon et al. (|1999|) 



an example, the fact that the (U, V) motion of the stream 
around IC 2451 A is very similar to the motion of an older 
(100-400 Myr) movi ng group inde pendent ly detected by 
Figuer as et al. ( 1997), A siain et al. ( 1999a ), Chereul et al. 
( |l999| ) and Sabas ( 1997 ) raises interesting questions. As 
a starting point, the vertical motions of these structures 
should enable us to confirm or reject any relationship. In 
this way, Comeron ( 1999 ) has recently reported a system- 
atic gradient in the vertical component of the velocity of 
those stars belonging to the Gould Belt along the galac- 
tic plane, which while being subtle is detectable in the 
Hipparcos astrometric data. 
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Fig. 16. Distribution of stars in the U-V plane - where a 
kernel estimator was used to indicate the lines of isocon- 
tours - for stars with 225 < Z < 285°, 100 < R < 300 pc 
and 30 < r < 60 Myr. The filled circles correspond to the 
position of the clusters in Table fj]. 



5. Conclusions 

A sample of O- and B-type stars with Hipparcos astro- 
metric data, radial velocities and Stromgrcn photometry 
- from which photometric distances and ages were com- 
puted - has been used to study the spatial distribution and 
the kinematics of the young star system in the solar neigh- 
bourhood. Several numerical simulations have allowed us 
to assess the robustness of our methods and to evaluate 
the biases induced by the observational constraints. 

The spatial distribution of the youngest and nearest 
group of stars is dominated by the presence of the Gould 
Belt. We found that this system extended up to 600 pc 
from the Sun and has an orientation with respect to the 
galactic plane of zq = 16-22° and £Iq = 275-295°, depend- 
ing on the distance and age intervals considered. For R < 
600 pc, roughly 60% of stars younger than 60 Myr belong 
to the Gould Belt. 

In the region with R > 600 pc, the stellar kinematics is 
dominated by the differential galactic rotation, since the 
Oort constants were found to be A = 13.0 ± 0.7 km s _1 
kpc" 1 , B = -12.1 ± 0.7 km s" 1 kpc" 1 , C = 0.5 ± 0.8 km 
s^ 1 kpc -1 and K — — 2.9±0.6 km s" 1 kpc~b In contrast, 
in the region with R < 600 pc, the Gould Belt dominates 
the kinematics of the youngest stars (r < 60 Myr), pro- 
ducing a decrease in the A and B Oort constants (A » 
6-8 km s" 1 kpc~\ B w -(21-14) km s" 1 kpc" 1 ) and an 
increase in C and K (C « 5-9 km s _1 kpc~ x , K k, 4-7 
km s _1 kpc -1 ). This peculiar kinematics was also found 



when those stars belonging to the Sco-Cen and Ori OBI 
complexes were eliminated. Therefore, these associations 
are not the only responsible for these peculiarities, a find- 
ing which seems to reinforce the suggestion by Guillout et 
al. ( |!998b| ) that the Gould Belt is a disk-like rather than 
a ring-like structure. 

A perfect agreement has been obtained when estimat- 
ing the age of the Gould Belt system from the spatial dis- 
tribution of the stars and from the study of the variations 
in the Oort constants with age. Taking into account the 
biases in the computation of individual photometric ages 
- stellar rotation and important uncertainties - we esti- 
mated an age of the Gould Belt inside the interval 30-60 
Myr. 

The study of the residual velocity field allowed us to 
estimate the cosmic dispersion for stars younger that « 
150 Myr. On the other hand, this residual velocity field for 
the youngest stars cannot be explained as an expansion 
from a point or a line. Moreover, the expansion motion 
classically attributed to the Gould Belt seems to be due 
to the nearest stars (i? Si 300 pc). In the region 300 < R < 
600 pc, we found that only Per OB2 has a clear residual 
motion away from the Sun. In the region with 100 Ss R Ss 
300 pc and 225 Si I Si 285°, two streams of stars with an 
age between 30 and 60 Myr have been found. One of these 
streams shares the motion of a Car (— HIP 45080) and IC 
2602, whereas the other follows the motion of NGC 2451 
A and Tr 10. 
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Appendix A: Simulations to check the structure 
analysis 

Simulations were performed to assess how incompleteness 
effects could change our conclusions on the Gould Belt's 
structure parameters. The critical questions to answer are: 

— Due to observational constraints older stars (r > 60 
Myr) have a small limiting distance (R « 400 pc). 
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If their spatial distribution shows an inclined struc- 
ture, as obtained for young and distant stars, could 
our method, which does not take into account incom- 
pleteness effects, be capable to detect it? 

— Related to the first point, for which scale height of the 
belts our method looses its statistical capability? 

— Are the available number of stars enough to under- 
take this study? Can simulations provide a realistic 
estimate of the errors in the derived structure param- 
eters? 

To answer these questions simulated samples were 
built by considering the following steps: 



From each real star we generated a pseudo-star with 
the same age, visual magnitude and projected distance 
(R cos b) . Its galactic longitude was randomly assigned 
and the distance to the galactic plane (z) simulated 
following an exponential distribution with scale height 
Z . 

Inside each age interval, the position of a fraction 
q = 0.50 of the generated pseudo-stars were rotated 
an angle «g = 20° around the Y axis (galactic rota- 
tion direction), that is adopting Qq — 270° for the 
ascending node of the Gould Belt. 
The process was repeated to generate samples with Z 
values ranging from 40 to 80 pc (Mihalas and Binney 
1981 quoted Zq = 60 pc for B-type stars). 



An example of the spatial distribution of the generated 
pseudo-stars is presented in Fig. Al , which can be com- 
pared with the distribution of the real stars in Fig. |l(| 
As seeked for, the same incompleteness effects are present 
in our simulated samples. The results after applying our 
resolution proces to the pseudo-stars with V < 7 and R 



< 600 pc are presented in Table Al, where in brackets we 
give the standard deviation for the 100 simulated samples. 

First, we confirm that the angular halfwidths (£g & n d 
£ g ) correctly reflect the growth of the scale height (Z ) of 
the simulated belts, the standard deviation of the differ- 
ents samples ranging from 2° to 5° . The q parameter is also 
well recovered (q = 0.50), though with a standard devia- 
tion as large as 0.13-0.17. On the other hand, although the 
values obtained for the angles defining the Gould Belt's 
orientation {iq, Qq) could indicate the presence of a small 
systematic trend when increasing Zq, probably due to the 



resolution process applied (Comeron 1992), it is always 
smaller than the standard deviation quoted. Looking at 
the interval 90 < r < 120 Myr, where the observational 
incompleteness is more accentuated, we realize that when 
the presence of the inclined structure (the Gould Belt) is 
present in the simulated sample, only below the 2a level 
we could obtain an iQ value as small as obtained from the 
real sample (ic = 3.4°). That is, there is a probability less 
than 5% to obtain inclinations about 4° in the worst case 
(Z - 80 pc). 



Table Al. Simulations on the Gould Belt's structural pa- 
rameters. Results obtained after averaging 100 simulated 
samples with input values: ?g = 20°, £!q = 270°, q = 0.50, 
with scale height: 40, 60 and 80 pc. Only pseudo-stars with 
R < 600 pc were considered. 



Z 






q 


cg n 


£g (°) 






T 


< 30 Myr 






40 
60 
80 


21.4 { 3.o) 
22.7(4.3) 

21.6( 5 . 2) 


268.0(8.7) 
268.5 { i3.a) 
268.5(20.3) 


0.50(o.i5) 
0.52 ( o.i7) 
0.51(o.i3) 


14.5(5.7) 
19.3(5.7) 

24.2(3.4) 


15-7(5.7) 
18.7(5.7) 

24.2(3.4) 






30 < 


t < 60 Myr 






40 
60 
80 


21-1(2.8) 
21-5(4.2) 

22.3(5.7) 


269.7(9.4) 

270.1 C i 5 . 9) 

268.1(31.6) 


0.50 ( o.i3) 
0.50(o.i4) 
0.49 ( o.i5) 


13.9(4.6) 
18.7(4.6) 

22.3(3.4) 


13.9(4.0) 

19.3(4.6) 

23.0(2.9) 






60 < 


r < 90 Myr 






40 
60 
80 


21-1(3.9) 

22.8( 6 .o) 
22.8(7.6) 


269.0(13.6) 
269.3 ( 22.5) 
262.7(50.7) 


0.50(o.i6) 
0.49 ( o.i8) 
0.50(o.i5) 


15.1(5.7) 
19.3(5.7) 

24.8(2.3) 


15-7(5.2) 

20.5(5.2) 
24.8(2.9) 






90 < 


t < 120 Myr 






40 
60 
80 


20.9(4.4) 

21-7(6.8) 

24.7(8.4) 


271.3(15.6) 

269.1(36.3) 

253.2(53.4) 


0.49 ( o.i7) 
0.48 ( o.i7) 
0.50(o.i6) 


15-1(5.2) 

20.5(5.2) 
25.5(2.3) 


16.3(5.2) 
21-1(5.2) 

25.5(2.9) 



Appendix B: Simulations to check the kinematic 
analysis 

Numerical simulations allow us to quantitatively evaluate 
the biases in the kinematic model parameters (Oort con- 
stants and solar motion components) induced by both, our 
observational constraints - irregular spatial distribution of 
the stars, incompleteness effects, availability of radial ve- 
locities, ... - and the presence of observational errors in 
the right hand side of Eqs. (§), (@) and (§), not consid- 
ered in our least square fit. Here we present the procedure 
followed to generate the simulated samples, the results de- 
rived when using them to solve Eqs. @, (f|) and (||), and 
finally, the quantification of the biases present in our real 
sample resolution. 

B.l. Process to generate the simulated samples 

To take into account the irregular spatial distribution of 
our stars and their actual observational errors, parameters 
describing the position of the each simulated pseudo-star 
were generated as follow: 

— From each real star we generated a pseudo-star that 
has the same nominal position (Rq, I, b) - not affected 
by errors - than the real one. 

— We assumed that the angular coordinates (l,b) have 
negligible observational errors. 

— If the star had a distance determination from the Hip- 
parcos parallax, the parallax error of the pseudo-star 
has a distribution law: 



e(Tr) 



(Bl) 
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Fig. Al. Example of the distribution of the simulated samples on the X-Z galactic plane. At the top, the pseudo-stars 
with an age less than 30 Myr; in the middle, those with an age between 30 Myr and 60 Myr; and at the bottom, those 
with an age larger than 60 Myr. Left: Gould and Galactic plane generated with Zq = 40 pc; right: both generated 
with Zq — 80 pc. 



where ov is the individual error in the parallax ttq 
(1/Rq) of the real star. From the value of n affected 
by the error, the simulated distance affected by error 
(R = 1/tt) was derived. On the other hand, if the star 
had a photometric distance determination, the error 
follows: 

_ i ( R-R \ 2 

e(R) = e *\ "b ) (B2) 

where is the individual error in the photometric 
distance of the star. 

To generate kinematic parameters we randomly as- 
signed to each pseudo-star a velocity (U, V, W) by assum- 
ing a cosmic dispersion (cry, oy, <Jw) and a Schwarzschild 
distribution: 

ip' v {U, V,W) = e n «u ) A 'v i A °w ) (B3) 



where (Uq,Vo,Wo) are the reflex of solar motion. These 
components were transformed into radial velocities and 
proper motions in galactic coordinates using the nomi- 
nal position of the pseudo-star (Ro,l,b). The systematic 
motion due to galactic rotation was added following Eqs. 
© ' @ an d (||) j obtaining the components (w ro , [i\ , /Zb ) 
for each star. Finally individual observational errors were 
introduced by using the error function: 



e(v T , m, (i h ) = 




where a Vi , er Ml and cr Atb are the observational errors of the 
real star. 

At the end of this process we had the following data 
for each pseudo-star: galactic coordinates (R, I, b), velocity 
parameters (v r , /j,b), errors in the velocity parameters 
(<7 Vl , <7 W , crp b ) and error in the trigonometric parallax (ov) 
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Fig. Bl. Fraction of pseudo-stars with radial velocity 
(top) and distribution of the stars with proper motion 
(blank histogram) and radial velocity (filled histogram) 
(bottom) as a function of proper motion for one of the 
simulated samples (3915 stars). Error bars were estimated 
from a Poissonian error distribution. 



or in the photometric distance (ctr). The simulated radial 
component of those pseudo-stars generated from a real 
star without radial velocity was not used, thus we imposed 
in the simulated sample the same deficiency in radial ve- 
locity data that is present in our real sample (Sect. |2.2f ). 
The fraction of pseudo-stars with radial velocity against 
proper motion is shown in Fig. Bl. This can be compared 
with Fig. H in Sect. |2.2| (the latter was made using the 
whole catalogue (6922 stars) whereas the former only con- 
tains the 3915 stars used in the fit of the Eqs. (||), (||) and 
(|J)). The systematic trend present in the real sample is 
well reproduced in the simulations. 



B.2. Results and discussion 

Following this scheme, two sets of 100 simulated samples 
were built, each one having the same number of stars than 
the real sample. The first set was built adopting K = 
km s _1 kpc -1 and used to derive the kinematic parame- 
ters in the distance interval 600 < R < 2000 pc. For the 
interval 100 < R < 600 pc we simulated t he e xpansion by 
imposing K = 5 km s -1 kpc -1 . In Table Bl we show the 



adopted kinematic parameters. For the cosmic dispersion 
we considered (a\j, cry, cw) = (8,8,5) km s _1 (see Sect. 
I- 



Case 1: A null error in distance was adopted. There- 
fore, the nominal distance of the pseudo-star is used 
(R = Rq) and only the effects of the errors on the 
radial velocity and proper motions are considered. No 
stars are rejected. 

Case 2: Errors in the radial velocity, proper motions 
and distance are considered. No stars are rejected. 
Case 3: As Case 2 but, to reproduce the real case, 
we rejected those stars with a residual velocity 3 times 
larger than the root mean square residual of the fit. 



Case 1 allows to study the effect of the errors in ra- 
dial velocity and proper motions, the incompletenes s o f 
our sample, the lack of radial velocity data (see Sect. 2.2 ) 
and the correlations between the different kinematic pa- 



rameters to be determined. As we can sec in Table Bl 



As can be seen in Table Bl , three different resolution 
processes were undertaken: 



no systematic bias is present, the difference between the 
adopted values and the obtained ones never exceeding 0.2 
km s _1 for the solar motion components and 0.2 km s -1 
kpc -1 for the Oort constants. 

When the error in the distance is considered (Case 2) 
the most noticeable effect is a clear bias in the A and B 
Oort constants. For A and the radial velocity resolution, 
a bias of +(0.6-0.7) km s _1 kpc -1 was obtained. On the 
contrary, from the proper motion equations this bias was 
— (0.7-0.8) km s™ 1 kpc -1 . For the combined resolution, 
and due to the larger number of proper motion equations 
(2iV against N), the bias was —(0.1-0.5) km s _1 kpc -1 , 
depending on the distance interval considered. For B con- 
stant, a bias of +(0.7-0.9) km s _1 kpc -1 was found from 
proper motion data. Apart from the distance errors, other 
effects contribute to the biases detected: the specific de- 
creasing distribution on distance of our real sample (well 
reproduced in the simulations) and the distance cut ap- 
plied in the resolution process (100, 600, 2000 pc). 

In Case 3, to mimic the real case, those stars with a 
residual velocity larger than 3 times the root mean square 
residual were rejected. In the simulated samples there are 
not high velocity stars, but in the reality it will be a few 
percentage of kinematically peculiar stars, or stars with 
non-well determined errors in the distance or in the veloc- 
ity components. We confirm that our rejection criterion do 
not induced any addicional bias, showing only an expected 
decrease in the x 2 statistic. 

Concluding, these simulations allow us to estimate 
which biases are expected in the kinematic parameters ob- 
tained from our real sample. As we can see, for 100 < R < 
600 pc our results for the combined solution can be biased 
on w —0.5 km s -1 kpc -1 in A Oort constant and w —0.8 
km s _1 kpc -1 in B, whereas for C and K is negligible. 
For the solar motion components a bias about 0.3-0.4 km 
s _1 can be present. For 600 < R < 2000 pc the bias for 
A, C and K is negligible while a positive bias of ~ 0.9 km 
s -1 kpc -1 is found for B. Again, for solar motion a bias 
of 0.3-0.4 km s _1 is expected in each component. 
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Table Bl. Mean Oort constants and residual solar motion for 100 simulated samples obtained soving Eq. (|3|) for 
radial velocities, Eqs. (|J) + (||) for proper motions and Eqs. (||) + + (j|) for the combined solution. The standard 
deviation for the 100 samples is shown in brackets. Units: A, B, C, K in km s _1 kpc -1 ; Uq, Vq, Wq, a in km s _1 . 
X 2 /JV eq is the value of x 2 divided by the number of equations. 



Radial velocities 







100 < R < 600 pc 






600 < R < 2000 pc 






Adopted 


Case 1 


Case 2 


Case 3 


Adopted 


Case 1 


Case 2 


Case 3 


A 


14.0 


13.9(i.8) 


14.7 (1 . 6) 


14.7(1.75 


14.0 


13.9 ( o.7) 


14.6 ( o.8) 


14.5 { o.8) 


C 


0.0 


0-2(i. 6 ) 


-0.1(1.7) 


-0.1(1.7) 


0.0 


—0.1(0.6) 


0-1(0.7) 


0-1(0.7) 


K 


5.0 


4.8(i.i) 


4.9 ( i.2) 


4.9(i.2) 


0.0 


—0.1(0.4) 


—0-3(0.9) 


—0-3(o.5) 


Uq 


9.0 


9.0(o. 5 ) 


9-1(0.4) 


9-1(0.5) 


9.0 


8-9(o.8) 


8-8(o.9) 


8-8(o.9) 


Vq 


12.0 


12-1(0.6) 


12-2(0.5) 


12-2(o.6) 


12.0 


12-1(0.7) 


11-6(0.7) 


11-6(0.7) 


Wq 


7.0 


7-2(i.o) 


7-2(i.o) 


7-3(i.o) 


7.0 


7.0 ( i.8) 


6-7(i.8) 


6-7(1.7) 


a 




8-4(0.2) 


8-5(0.2) 


8-4(o.3) 




8-3(o.3) 


9-2(0.4) 


9-0(0.4) 






1.00 


1.00 


0.97 




1.00 


1.11 


1.06 


Proper motions 






100 < R < 600 pc 






600 < R < 2000 pc 






Adopted 


Case 1 


Case 2 


Case 3 


Adopted 


Case 1 


Case 2 


Case 3 


A 


14.0 


14.2 ( o.8) 


13.2(o.8) 


13.2(o.8) 


14.0 


14.0 { o.5) 


13.3(o.5) 


13.3(o.5) 


B 


-12.0 


— 12.0(o.7) 


— 11-3(0.7) 


-11-3(0.7) 


-12.0 


— 12.0(o.4) 


-11-1(0.4) 


— 11-1(0.4) 


C 


0.0 


0-1(0.8) 


-0.1(0.8) 


—0.1(0.8) 


0.0 


0.0(0.5) 


— 0-1(0.6) 


—0-1(0.6) 


K 


5.0 


5-2(i.8) 


4-8(i.8) 


4-8(i.9) 


0.0 


0.1(2.1) 


—0.2(2.1) 


—0.3(2.1) 


Uq 


9.0 


9-0(0.2) 


8-6(o.2) 


8-6(0.2) 


9.0 


9-0(0.6) 


8-7(o.7) 


8-7(o.7) 


Vq 


12.0 


12-0(o.2) 


11-4(0.2) 


11-4(0.2) 


12.0 


12-1(0.6) 


11-4(0.7) 


11-4(0.7) 


Wq 


7.0 


7.0 ( o.i) 


6.6(o.i) 


6.6(o.i) 


7.0 


7-1(0.3) 


6-7(0.3) 


6-7(0.3) 






6-4(o.i) 


6.6(o.i) 


6.6(o.i) 




7-6(o.2) 


7-8(0.2) 


7.6(0.2) 


X 2 /N cq 




1.00 


0.94 


0.91 




1.00 


0.91 


0.88 


(X7^e q ), 




1.00 


0.93 


0.90 




1.00 


0.87 


0.83 


(X 2 /^ q )b 




1.00 


0.94 


0.91 




1.00 


0.96 


0.93 


Combined solution 






100 < R < 600 pc 






600 < R < 2000 pc 






Adopted 


Case 1 


Case 2 


Case 3 


Adopted 


Case 1 


Case 2 


Case 3 


A 


14.0 


14.1(0.7) 


13.5(o.7) 


13.5(o.7) 


14.0 


14.0 ( o.6) 


13.9(o.5) 


13.9 { o.5) 


B 


-12.0 


— 12.0(0.7) 


— 11-2(0.7) 


— 11-2(0.7) 


-12.0 


— 12.0(0.6) 


-11-1(0.4) 


— 11-2(0.4) 


C 


0.0 


0-1(0.7) 


0-2(0.7) 


0-2(0.7) 


0.0 


-0.0(0.5) 


0-1(0.4) 


0-1(0.4) 


K 


5.0 


5-0(0.8) 


5-1(0.8) 


5-1(0.9) 


0.0 


-0.1(0.5) 


— 0.2(0.4) 


—0.2(0.4) 


Uq 


9.0 


9-0(0.2) 


8-7(0.2) 


8-7(o.2) 


9.0 


9.0(0.7) 


8-8(o.5) 


8-7(o.5) 


Vq 


12.0 


12-0(0.2) 


11-6(0.2) 


11-6(0.2) 


12.0 


12-1(0.6) 


11-6(0.5) 


11-6(0.5) 


Wq 


7.0 


7-0(0.1) 


6-7(o.i) 


6-7(o.i) 


7.0 


7-1(0.4) 


6-7(0.3) 


6-7(0.3) 


a 




6.7(o.i) 


6-7(o.i) 


6.6(o.i) 




10-2(0.2) 


8-1(0.2) 


8.0(o.2) 






1.00 


0.95 


0.92 




1.00 


0.97 


0.93 


(x 2 /^ cq ) r 




1.00 


1.00 


0.97 




1.00 


1.12 


1.08 


(x 2 M q )i 




1.00 


0.93 


0.90 




1.00 


0.98 


0.84 


(x 2 /^v cq )b 




1.00 


0.94 


0.91 




1.00 


0.96 


0.93 



